arXiv:1504.01224vl [cond-mat.soft] 6 Apr 2015 


Condensed Matter Physics, 2015, Vol. 18, No 1, 13606:1 -[s] 
DOI: 10.5488/CMP.18.13606 




http://www.icmp.lviv.ua/journal 


Heat capacity of liquids: A hydrodynamic approach 

T. BryPl, T. ScopigncP, G. RuocccP 


^ Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, 

1 Svientsitskii St., 79011 Lviv, Ukraine 

^ Institute of Applied Mathematics and Fundamental Sciences, Lviv Polytechnic National University, 
79013 Lviv, Ukraine 

^ Dipartimento di Fisica, Universita di Roma La Sapienza, 1-00185, Roma, Italy 
^ Center for Life Nano Science OSapienza, Istituto Italiano di Tecnologia, 

295 Viale Regina Elena, 1-00161, Roma, Italy 

Received March 3, 2015 


We study autocorrelation functions of energy, heat and entropy densities obtained by molecular dynamics sim¬ 
ulations of supercritical Ar and compare them with the predictions of the hydrodynamic theory. It is shown 
that the predicted by the hydrodynamic theory single-exponential shape of the entropy density autocorrelation 
functions is perfectly reproduced for small wave numbers by the molecular dynamics simulations and permits 
the calculation of the wavenumber-dependent specific heat at constant pressure. The estimated wavenumber- 
dependent specific heats at constant volume and pressure, Cvik) and Cp(fc), are shown to be in the long- 
wavelength limit in good agreement with the macroscopic experimental values of Cv and Cp for the studied 
thermodynamic points of supercritical Ar. 
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1. Introduction 

Thermodynamics of the liquid state of matter was successfully explained on the basis of thermody¬ 
namic perturbation theory (TDPT) more than 40 years ago (im). Together with the integral equation 
theory of static structure of liquids and computer simulations, it is the basis of current exploration of 
physics and chemistry of real liquid systems. Modern versions of TDPT (5][6j were recently developed on 
the basis of treatment of the interaction interatomic potential as a sum of an improved reference system 
(hard-sphere- or soft-sphere-like repulsion plus short-range Yukawa-hke attraction) and a perturbation 
due to the remaining tail of the interparticle interaction. 

It is also possible to represent the thermodynamics of the liquid state via time-dependent processes 
because dynamic processes in condensed matter systems are tightly connected with their thermodynam¬ 
ics. However, in contrast to solids, where thermodynamics can be easily defined by contributions from 
static lattice and small vibrations around the lattice sites in terms of phonon excitations (7l, in fluids 
such an approach is not applicable due to the dynamic disorder and the absence of stable local energy 
minima for atomic particles. Furthermore, in liquid state there exist various relaxation processes which 
are absent in the solid state. As an example one can consider a relaxation connected with diffusivity of 
local temperature. In hquid system without any external temperature gradient, this relaxation process 
is caused by local temperature gradients which emerge due to adiabatic propagation of the most long- 
wavelength acoustic excitations. The acoustic modes, which propagate with smaU wavenumbers due to 
local conservation laws, in fact create adiabaticaUy compressed and decompressed local regions in liquid 
with different local temperatures. Such a complex origin of collective dynamics in fluids which is com¬ 
pletely different from the elastic sound propagation in solids is tightly connected with the phenomenon 
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of viscoelasticity of liquids (8). Furthermore, the acoustic excitations in liquid state are only of longitu¬ 
dinal polarization (transverse long-wavelength acoustic modes are not supported in liquids, there can 
exist only short-wavelength shear excitations with a propagation gap in the long-wavelength region) and 
have non-zero damping. The dispersion law of longitudinal and transverse collective excitations in liq¬ 
uids cannot be considered as the linear in wave numbers (9 ][To 1 as this is taken in the Debye model of 
the heat capacity of solids. All these facts must be reflected in the construction of liquid thermodynamics 
represented via the dynamic processes. 

Heat capacity being one of the important quantities in thermodynamics of the liquid state of matter 
can be easily calculated from the temperature fluctuations in the liquid system (TJ}. This is by date the 
most simple way of estimation of the specific heat at constant volume Cv for realistic liquids from com¬ 
puter simulations. Schofield (T2|[T3 has shown how the generalized wavenumber-dependent thermody¬ 
namic quantities can be estimated from static correlators of different statistically Independent dynamic 
variables m. The generalized wavenumber-dependent thermodynamic quantities were studied later on 
for pure irsHTTl and binary liquids in classical and ab initio (T9l molecular dynamics simulations. 

Perhaps the first application of the hydrodynamic theory to collective dynamics of liquids was re¬ 
ported in the paper by Landau and Placzek ^0), while systematic theoretical studies of liquid dynamics 
with hydrodynamic equations have started since the seminal studies by Mountain (21). In 1971 Cohen 
et al. derived analytical expressions for hydrodynamic time correlation functions for pure and binary 
liquids (22), which contained the so-called “non-Lorentzian corrections” and satisfied all the sum rules 
existing within the hydrodynamic treatment (23). Molecular dynamics (MD) simulations being a perfect 
tool for exploration of time-dependent correlations in realistic liquids and solids enable one to check the 
predictions of the hydrodynamic theory. 

In this study we report the behaviour of autocorrelation functions of energy, heat and entropy densi¬ 
ties obtained by molecular dynamics simulations of supercritical Ar at two densities. By date there were 
no MD studies of the entropy density autocorrelation functions of pure liquids, for which the hydrody¬ 
namics predicts purely single-exponential decay with time in the long-wavelength region. Besides, the 
zero-time value of the entropy density and heat density autocorrelation functions should be connected 
with the wavenumber-dependent specific heats Cp and Cy, respectively. We will perform a check of the 
long-wavelength limits of the Cp(A;) and Cv(fc) with the macroscopic experimental values for these ther¬ 
modynamic quantities l24l . The remaining part of the paper is organised as follows. In the next section 
we give the predictions of the hydrodynamic theory on the studied time correlation functions. In section 3 
we supply details of the MD simulations, and section 4 reports obtained the results on the autocorrelation 
functions of energy, heat and entropy densities obtained by molecular dynamics simulations of supercriti¬ 
cal Ar, their corresponding time-Fourier transforms, and the estimated dependencies for the generalized 
specific heats Cp(fc) and Cv(fc). The last section contains conclusions of this study. 


2. (Hydrodynamic approach 

Hydrodynamic approach to collective dynamics of liquids is valid only on macroscopic length and 
time scales where the most long-time relaxation processes in continuum are considered. The macro¬ 
scopic equations in the hydrodynamic theory are in fact the local conservation laws for the densities of 
particles, of total momentum and of energy. For the case of longitudinal dynamics, the only relevant on 
macroscopic scale contributions to the hydrodynamic time correlation functions (8) are the ones coming 
from acoustic longitudinal propagating modes with the linear in wave numbers k dispersion law 


Whyd(fc) = Csfc 


and quadratic in k damping 

O-hydffc) = r/c^. 

Here Cs is the adiabatic speed of sound and the damping coefficient F = (Dl + {j- 1)Dt) /2 depends on the 
longitudinal kinematic viscosity Dl, the thermal dlffusivity Dy and the ratio of specific heats j - Cp/Cy. 
The only relevant relaxing mode for the longitudinal dynamics on macroscopic scales Is the one coming 
from thermal relaxation, i.e., connected with the diffusivity of local temperature Dy. 
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The analytical expressions for the hydrodynamic time correlation functions, which we will study by 
MD simulations, have in general a three term form: an exponential function connected with thermal 
relaxing mode and two oscdlating terms coming from the acoustic excitations (^. The heat density au¬ 
tocorrelation function reads 


Fhh(k,t]^ Ahhe + [BhhC0SCskt + Dhhik]sincskt]e , ( 2 . 1 ) 

where the amplitudes of mode contributions are constrained by the lowest-order sum rule 

Ahhik) + Bhhik)^kTiT^Cy[k). ( 2 . 2 ) 

Note that the mode contributions A}jh and B^h are constants in the long-wavelength limit, whUe the am- 
plitide of the “non-Lorentzian” term D^hik) linearly depends on k. However, outside the hydrodynamic 
region when one applies MD simulations the amplitudes A^f, and B^f, become fc-dependent making pos¬ 
sible estimation of the wave-number dependent specific heat. Similar Is the expression for the energy 
density autocorrelation function: 

Fge{k,t) - Aeee~^'^^ ^ + [BggCosCskt + Deeik)smc^kt]e~^^ \ (2.3) 

while for the entropy density autocorrelation function, the hydrodynamics predicts a single exponential 
form for pure liquids: 

Fss{k,t)^Asse-^^^"\ (2.4) 

with 

Assik) ^ k^Cpik) . (2.S) 

Hence, the wavenumber-dependent specific heats Cv(fc) and Cp{k) are connected with the static heat- 
density and entropy-density autocorrelation functions, which in fact are the zeroth moments of the cor¬ 
responding spectral functions Shhik,co) and Sss{k,cD). In the long-wavelength limit, the wavenumber de¬ 
pendent specific heats tend to their macroscopic values Cv and Cp. The advantage of the hydrodynamic 
approach lies in a possibility to estimate the contributions from relaxing and propagating modes to cor¬ 
responding thermodynamic quantities f25l . 


3. Molecular dynamics simulations 

We performed molecular dynamics simulations for supercritical Ar at temperature T - 280 K and two 
densities 750 kg/m^ and 1464.5 kg/m^ using a system of 2000 particles Interacting via the Woon potential 
(26] with the parameters taken from (23. The cut-off radius for the effective two-body potential was 12 A. 
This potential permits a nice reproduction of the thermodynamic and dynamic quantities of supercritical 
Ar by MD simulations as it was shown in (28]. 

The simulations were performed in the microcanonical ensemble with the time step of 2 fs. The du¬ 
ration of production runs was of 360 000 time steps. Every sixth configuration was used for the sampling 
of dynamic variables of particle density, momentum density, energy density and first time derivative of 
the momentum density. The averages of static and time correlation functions over aU possible directions 
of different wave vectors with the same wave number were performed. 

During the MD simulations we sampled the spatial Fourier components of hydrodynamic variables 
of: particle density 

1 ^ 

n{k, t) = —= y , (3.1) 

VNti 

density of longitudinal momentum 


J^ik, t] = 


N 


k\fNi^i 




„-ikri(t) 


(3.2) 
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and energy density 


1 ^ 

e{k, t) = ——YetWe 


ikr,(f) 


(3.3) 


where r;(f), v,(f) and Ei{t] are the trajectories, velocities and single-particle energies of the i-th particle 
Ql- The single-particle energies were calculated as follows: 


Eiit) = 






where <l)(r) is the effective pair potential used in the simulations. 

The above dynamic variables were used to calculate the hydrodynamic time correlation functions as 
well as the heat density autocorrelation functions 


Fhh^k, t) - {h[k, t]h* [k, t - 0)) 


(3.4) 


with 

fne 

h{k, t) - e{k, t) - —— n{k, t), (3.5) 

Inn 

where the brackets mean the ensemble average and static averages read as follows fneik) - {n{k)e{-k)) 
and fnnik) = {n[k)n{-k)). 

In order to study the entropy density autocorrelation functions Fssik, t), we sampled the density of 
the longitudinal component of stress tensor Yik, t) via QH 

Y{k, t) = -ika^ik, t), 


where the overdot denotes the time derivative of the corresponding dynamic variable. Hence, the dy¬ 
namic variable of entropy density can be sampled from MD simulations as Ill4| 


s[k, f) = 


1 

f 


e(k, t) 


^n[k,t) 


with the corresponding static averages fg„L{k) - (e{k)a^(-k)) and fuL„{k) - (cr^(A:)n(-fc)). 


4. Results and discussion 

Energy density autocorrelation functions Fgelk, t) - (e{k, t)e* (fc, t - 0)) calculated in MD simulations 
are shown for three wavenumbers in small-fc region in figure The time dependence of these time cor¬ 
relation functions is very similar to the relaxation of the heat density autocorrelations shown in figurej^ 




Figure 1. (Color online) Energy density autocorrelation functions Fee (k, t] at three wave numbers as 
calculated from MD simulations for supercritical Ar at T = 280 K and densities 750 kg/m^ (a) and 
1464.5 kg/m^ (b). The time scales t are equal to 3.701 ps (a) and 2.961 ps (b). 


13606-4 





















Heat capacity of liquids 


A typical feature is an increasing antplitude and a decreasing frequency of the damped oscillations in 
the shape of Fgeik, t) and t) towards smaller wave numbers, which is typical of the contribution 

coming from acoustic excitations. The relative strength of the relaxing and oscillating contributions in 
Fhh^k, t) was estimated in the long-wavelength limit in as the inverse of ( 7 - 1 ), i.e., the inverse of 
the Landau-Placzek ratio known for the density-density time correlation functions. Hence, for the case of 
the two thermodynamic points of the supercritical Ar studied here and having the experimental ratio of 
specific heats 7 of 2.30 and l.SS (24) for the low- and high-denslty studied systems, one should observe 
much smaller relaxing contribution for the low-density state. Indeed, in figure (a) the oscillations in 
the shape of Fhh{k, t) even make the heat density autocorrelation function negative close to wt ~ 1 for 
the smallest possible in these simulations wave number, while in figure [^(b), due to the much stronger 
contribution from the relaxing thermal mode, the F^f,{k, t) is positive for all times. It is important that 
the dimensionless functions F^hik, t] have the zero-time values tending to the macroscopic values of Cy 
due to l |2.2| . Indeed, in figure[^(a) the zero-time values tend to the macroscopic value of Cy = 1.775 1^ . 
while for the hlgh-density state they tend to the value Cy = 2.34. 





Figure 2. (Color online) Heat density autocorrelation functions F^h^k, t) at three wave numbers as calcu¬ 
lated from MD simulations for supercritical Ar at T = 280 K and densities 750 kg/m^ (a) and 1464.5 kg/m^ 
(b). The time scales t are equal to 3.701 ps (a) and 2.961 ps (b). 




Figure 3. (Color online) Entropy density autocorrelation functions ^SS (k, t) at three wave numbers as 
calculated from MD simulations for supercritical Ar at T = 280 K and densities 750 kg/m^ (a) and 
1464.5 kg/m^ (b). The time scales t are equal to 3.701 ps (a) and 2.961 ps (b). 

The entropy density autocorrelation functions Pss [k, t) calculated In MD simulations for two den¬ 
sities of the supercritical Ar (figure show for the smallest wave numbers a simple relaxation shape 
without any oscillating contribution. This is in complete agreement with the prediction of the hydrody¬ 
namic theory |2.4|. The zero-time values of F^^ik, t) (2.5| should tend in the long-wavelength limit to the 
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macroscopic value of Cp. Indeed, one observes in figures[^(a) and (b) that in agreement with this hydro- 
dynamic prediction, the zero-time values of the normalized functions tend to the macroscopic values of 
Cp reported in the NIST database: 4.097 for the low-density state and 3.62 for the high-density state. 

The time-Fourier transformed energy density (e-e), heat energy density (h-h) and entropy density (s-s) 
autocorrelation functions, Sa [k, oj), i- e, h, s are shown for the two studied densities of the supercritical 
Ar in figure]^ The See{k,cD] and Shhik,(jj] have the typical of the dynamic structure factors shape with 
the central peak due to thermal relaxation and the side peak due to propagating acoustic excitations. 
The ratio of the integral intensities of the central and two side peaks known as the Landau-Placzek- 
type ratio was obtained in (23 for the case of Shh{k,(D) and is equal to (y- 1 )“^. The side peak of the 
See(k,(jj] and Shh{k,(D) is the consequence of the corresponding oscillating contributions to the Feeik, t] 
and Ffiijlk, t), respectively. It is remarkable that the calculated spectral function Sssik.cj) shows only the 
one-peak shape with the central peak due to thermal relaxation (line connected star symbols in figure]^. 



Figure 4. (Color online) Spectral functions energy-energy Seeik,(i)), heat-heat ShhikyCx)), entropy-entropy 
Sssik,o)) at the smallest available in MD simulations wave numbers for supercritical Ar at T - 280 Kand 
densities 750 kg/m^ (a) and 1464.5 kg/m^ (b). The time scales t are equal to 3.701 ps (a) and 2.961 ps (b). 
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Figure 5. (Color online) Generalized wavenumber-dependent specific heats at constant volume Cv(fc) and 
constant pressure Cp(fc) for supercritical Ar at T = 280 K ad densities 750 kg/m^ (a) and 1464.5 kg/m^ (b). 

The asterisks at fc = 0 correspond to macroscopic values of Cv calculated from the temperature fluctua¬ 
tions im . 

The zero-time values of the time correlation functions F^hik, t] and Fgsik, t) permit estimation of the 
wavenumber-dependent specific heats Cv(fc) and Cp(A:), which in the long-wavelength limit tend to their 
macroscopic values. In figurej^we show the wavenumber dependencies of Cyik) and Cp(fc) for the two 
densities of the supercritical Ar at temperature 280 K. The asterisks at fc = 0 show the values of macro- 
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scopic values of Cy obtained from the standard expression via the thermal fluctuations in the fluid system 
El. One can see that the reported wave-number dependencies tend quite well to the experimental val¬ 
ues for the supercritical Ar taken from the NIST database (24l: Cv=1.77SfcB and Cp=4.097A;B for the density 
of 750 kg/m^ and Cv=2.34A:b and Cp=3.62fcB for the density of 1464.5 kg/m^. 

5. Conclusions 

In the MD simulations we calculated the autocorrelation functions of energy density, heat density 
and energy density and compared them for the smallest available wave numbers with the predictions 
of the hydrodynamic theory. Our results give evidence of a nice correspondence of the hydrodynamic 
description of time correlation functions and their zero-time values, which made it possible to estimate 
the wavenumber-dependent specific heats Cv(fc) and Cp(A:), with the MD-based estimation of their macro¬ 
scopic values and experimental values for specific heats of supercritical Ar taken from the NIST database. 

These results are very important from the viewpoint of analysing the contributions from thermal 
relaxation and collective excitations to the specific heats of liquids. Our results and the hydrodynamic 
approach make evidence of the definitive role of the thermal relaxation in calculation of the specific heats 
of fluids when it is estimated via dynamic processes. These results will urge a further understanding 
of the dynamics of the supercritical state of matter, which shows many interesting features observed 
recently in the scattering experiments and MD simulations on supercritical Argon (29] |30| and Oxygen 
[3ll. 
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Ten/iocMHiCTb piAi/iH: riApoAi/iHaMiHHi/iii niAxiA 

T. Bpi/ikl® T. CKoniHbMcP, /],x. PyoKKcP^I 

^ iHCTHTyr (fiisuKM KOHfleHCOBaHnx CMCreM HAH yKpaiHi/i, 79011 /IbBie, VKpaTHa 

^ iHCTi/iTyr npMK./iaflHoT MareMaTHKi/i ra c|)yHflaMeHTa/ibHMx Hayx, HaL(ioHa./ibHMM YHiBepcuTeT "/IbBiBCbxa 
no/iirexHiKa", 79013 /IbBiB, yxpaiHa 

^ (Piai/itHi/iM (()aKy./ibTeT, PuMCbKMw yHiBepcMTeT"La Sapienza", 1-00185, Pi/im, lTa./iia 

^ peHTp 6io-HaHO Hayx Sapienza, lTa./imcbKnPi TexHO.fioriHHni1 iHCTHTyr, 295 syn. Kopo/ieen E.fieHi/i, 1-00161, 

Pmm, Ira/iia 

Mm AOcaiflxycMO aBTOKopeatipiCiHi cfjyHKpii rycTMH eHeprii', Ten.na ra eHTponii, oOnucaeHMX npn MOfle/iiOBaH- 
Hi MeroflOM M0./ieKy./itipH0T flUHaMiKi/i pfla HaAKpniMtHoro At ra nopiBHKDCMO Tx is nepefleaneHHtiMM riflpofli/i- 
HaMiHHOi reopiT. noxasaHO, tpo nepeflOanena riflpofli/iHaMinHOio reopieio oflHoexcnoHeHTHa c|)opMa bbtoxo- 
pe.fiaL(iMH0T cfiyHxpii' rycn/iHi/i eHiponii' nyflOBO BiflTBopioeTbCB npn Ma./ii/ix XBi/i/ibOBHx HMC/iax MO/iexy.fiapHOK) 
fli/iHaMixoK) ra floSBO/ine poapaxynox 3a./iexHOi Bifl xBM.fibOBoro HMC.fia omtomoi Ten/ioewHOCTi npn CTa.fiOMy tm- 
cxy. OrpnMani nnroMi Ten.noeMHOCTi, 3a./ie>XHi Bifl XBM./ibOBOro nucaa, npi/i CTa.noMy o6'eMi ra n/icxy, Cp(k] and 
Cp[k), y flOBtoxBM/ibOBiPi rpaHupi flo6pe y3rofl>xyK)TbCti 3 MaxpocxoninHMMU excnepMMeHTa.fibHMMM ananeH- 
HtiMM Cv and Cp fl/ia flOC.fiifl>xeHMX repMOflUHaMinHnx ronox HaflxpMTMHHOro At. 

KnioMOBi cnoBa: repMOflUHaMiKa, ren/ioeMHicTb, ri^po^i/iHaMi^Ha reopifi 
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